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REMARKS/ARGUMENTS 

In response to paragraph 2 of the Office Action (page 2), applicants have now removed all the new 
matters (equations, functions, matrices, etc) as the Examiner required. Applicants submit that the new 
Specification includes no new matter, and respectfully request withdrawal of the objection. 

In response to paragraph 3 of the Office Action (page 3), 

In response to paragraph 4 of the Office Action (page 3), applicants haye added Figure 1 and Figure 2, 
so as to make the specification easier to understand, please note that 

(1) Figure 1 and Figure 2 only show well-known concepts in automatic control; and 

(2) Although the present invention can be clearly described in words as: 

"This invention chooses tlie best Values for the tuning parameters in a PID controller 
or a linear controller in such i^AA/ay that the largest absolute value of all poles of said 
discrete-tima closed-lodfrtra^ said set-point r(k) to said process 

variable v(k) is minimized subject to. if any, user-specified constraints on one or more 
of the tuning pararrteters." (see the Specification or the Claims) 

this invention cannot be described by drawings, because it contains abstract mathematical concepts 
(see the above underlined words) that cannot be described by drawings. However, anyone skilled in the 
art can easily understand the invention as described in the Specification and Claims and can easily 
carry out this invention by using well-developed optimization algorithms in mathematics or 
well-developed computer programs such as the "Optimization Toolbox" in Matlab developed by The 
MathWorks Inc. (Matlab is a registered trademark of The MathWorks Inc.). The applicants do not need 
to teach these well-known prior arts. 
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In response to paragraph 5 of the Office Action (page 4), applicants submit that the new specification is 
now in full compliance with the layout suggested by the Examiner. 

In response to paragraph 6 of the Office Action (page 8), applicants submit that the new dependant 
claims now refer to other claims in the alternative only and comply with 37 CFR 1.75 (c). Applicants 
respectfully request re-examination of the new claims and withdraw of the objections. 

In response to paragraphs 7, 8 and 9 of the Office Action (pag£s 8-9), applicants have removed the 
equations, transfer functions, and matrices as required by the Examiner. Anyone skilled in the art 
knows how to find the closed-loop or open-loop transfer functions, the pbles of a transfer function, the 
definition of PID controllers and linear controllers and their structures, and the use of well-developed 
optimization algorithms in mathematics and optimization corirputer programs such as the "Optimization 
Toolbox" in Matlab developed by The MathWoi'ks lnc> to solve the optimization problem formulated in 
the present invention (Matlab is a registered trademark ;^^;^iyiathWorks Inc.). The applicants do not 
need to teach these well-knowri prior arts . For example, if an invention needs the calculation of the 
value of sin(x), the inventor does ribfe;need^tp teach the- definition of sin(x), neither does he need to 
teach the various mathematical rriethods, algorithms or computer programs for finding the value of 
sin(x), since these are all, well-known prior arts/ In Matlab, the function 'Tminimax" can be used to 
directly solve the optimizatibni problem i(minimax problem) formulated in this invention (see the newly 
cited book "Optimization Toolbox User's Guide" authored by Coleman et al. that comes with the 
commercial software product "Optimization Toolbox" developed by The MathWorks Inc). Matlab and the 
associated "Optimization Toolbox" are well known among people skilled in the art. Therefore, applicants 
do not need to teach how to minimize the maximum of absolute values of all poles of a discrete time 
closed-loop transfer function from a set point to a process variable. 

Therefore, applicants respectfully submit that the new claims have overcome the defects as pointed out 
by the Examiner and are in full compliance with 35 U.S.C. 112 since anyone skilled in the art can 
implement this invention without any difficulty. Applicants respectfully request re-examination of the 
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new claims and withdrawal of the rejections and criticisms as stated in paragraphs 7, 8 and 9 of the 
Office Action. 

In response to paragraph 10 of the Office Action (page 9), applicants revised the specification so that it 
now relies on the prior art of record cited. 

Applicants have read through the following new reference cited by the Examiner: 

A. Cichocki et al., "Neural Networks for Solving Systems; of Lihear Equations II: Minimax and 
Least Absolute Value Problems", IEEE Transaction on Circuits and Systems II: Analog and 
Digital Signal Processing, Volume 39, Issue 9, Sept. 1992, Pages 619-633. 

A careful study of this prior art shows it deals with the problem of solving linear equations by minimizing 
the largest absolute error using the neural networks ^qhnology. It^is not apparent that this method can 
be used to solve the minimax problern formulated in this itiyention since the minimax problem in this 
invention is totally different, Hovyever, there are many wefl-developed optimization methods in 
mathematics that can directly solv£\ihe nnrinirnax probieto formulated in this invention (see, e.g., the 
cited book "Minimax and Applications" authored by ; Du et al. and published by Kluwer Academic 
Publishers, Dordrecht, in 1995). Indeed many ^computer program packages such as the "Optimization 
Toolbox" in Matlab developediby The Math Works Inc. can solve the minimax problem formulated in this 
invention without any difficulty. Matlabraild its various toolboxes such as the "Optimization Toolbox" are 
very popular among people working in automatic control area. All these well-developed optimization 
methods are well-known prior art and do not need to be taught in this invention. 

Applicants respectfully submit that the new specification now relies on the cited references and 
therefore complies with the requirements in paragraph 10 of the Office Action. 

Conclusion 



PACE 10/32 * RCVD AT 7/12/2004 1:26:18 AM [Eastern Daylight Time] * SVR:USPTO-EFXRF-1/0 * DNIS: 8729306 * CSID: 7804245074 • DURATION (mm-ss): 17-44 



07/11/2004 23:25 



Dan Liu 



PAGE 11/32 



Application/Control Number 10/069,364, Art Unit 2121 Page 10 

In order to comply with the Office Actions, Applicants have re-written the specification and claims. For 
all the aforesaid reasons applicants submit that the new specification and claims do not contain new 
matters. Drawings are added. The layout of the specification is now as suggested by the Examiner. The 
specification and claims enable anyone skilled in the art to carry out the invention without any difficulty. 
Applicants respectfully submit that this application is now in full condition for allowance. 

However applicants are willing to accept any suggestions that the Examiner believes necessary in 
order to avoid further proceedings. 

Very respectfully, 



Danyang Liu Wuyang Liii 

#1 08, 1 061 5-111 Street P.O. Box 1 1 2 : 

Edmonton, Alberta T5H 3G1 , % 316, San Jihg Lu F 

Canada Napchang, P. F?, China 
Telephone: (780) 424-5074 
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VERSION WITH MARKINGS TO SHOW CHANGES MADE 
In the Specification: 

Please cancel the old specification in record and replace it with the new Specification (underlined) as 
follows: 

Doscript i on 

T I TLE OF THE I NVENT I ON 

Method s of D e s i gning Opt i ma l P I D Contro ll ors 



TECHN I CAL F I ELD OF THE I NVENT ION 




Th i c invont i on ro l atoc to tho dosign of tho structure of a multiyariablo P I D cbntro ll or and tho optima l 
choico of i ts PID paramotors. 

BACKGROUND OF THE I NVENTION 

A trad i t i ona l PID contro ll or is usod ta%ntro^ Tho procoss variab l o (PV) goos i nto 

tho PLD control l or, which calculatb&tho cbntrollor output (CO) according to a P I D control oquat i on. This 
CO is then convortod tff gn ana l og signal, which is sent to tho proc e ss so that tho said PV can track a 
usor spoc i fiod va l uo ca ll od sot point (SP),Tho sa i d SP can chango w i th t i mo. Tho porformanco of a P I D 
contro ll or doponds on tho choicb bf ite throo PID paramotors. For i ndopondont form of P I D contro ll ors 
th e s e throo P I D paramotors aro th e proport i onal ga i n Kp, the i nt e gr al gain K i , and th e d e r i vat i v e ga i n 
Kd. For d e p e nd e nt form of P I D control le rs thoso throo P I D paramotors are tho gain K, i ntegra l time T i T 
and dorivat i vo t i mo Td. I n trad i t i ona l PID contro ll or s tho s a i d PV, SP, CO, and PED paramotors aro a ll 
scalars. Wo cal l this k i nd of PID contro ll ors tho s i ng l o -i nput singlo - output (S I SQ) P I D contro ll ors. Tho 
Z i og l or - N i ohols P I D contro llo r tun i ng m e thod i s th e m a jor on e of th e many m e thods for f i nd i ng th o 
va l ues of P I D parameters. 

DETA I LED DESCR I PT I ON OF THE INVENTION 
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I n this i nvont i on tho S I SQ PID controHor i s oxtondod to tho mult i plo - input mu l t i p l o - output (MIMO) P I D 
contro ll or that has n procos s var i ab l es PV1, PV2, = = = , and PVn and m control l or outputs CQ1 , 
CQ2, . . . , and COm, whoro m and n aro po si t i vo i ntogors, Correspond i ng to PV1 , PV2, .... and PVn 

th e r e ar e n s e t points SP1 , SP2, . . . , and SPn. I n this cas e PV bocomos a vector w i th PV1, PV2 

and PVn be i ng i ts f i rst, second, .... and n - th componont, CO becomes a vector with CQ1 , CQ2, .... 
and COm bo i ng i ts first, soGond, . . . , and m - th componont, SP bocomos a voctor with SP I , SP2, 
and SPn bo i ng i ts fir s t, s econd, . . . , and n - th componont, and tho P I D contro l oquation bocomos 
CO(k) - CO(k - 1)+KrSP(k)*T+Kra(k,1)+K2 ,,f a(k,2) +. . . +Kj*a(k,j), wh o r e k i s th e d i screte t i me, T i s tho 
campling por i od, j is a pos i tivo intogor, K1 , K2, = . . , Kj aro m by n P I D paramotor s , a(k,1)-[ - PV(k)]*T, 
and a(kj)-[a(k,j - 1) - a(k - 1 j - 1)]/Tfor j> or~2. I t i s i mportant tp fi6to that v., 

1. an M I MO PID control l or i s ab l e to tak e i nto account "the i nteract i on amorigVthe n proc e ss variab le s 
and m control le r outputs, which can not bo achipvod by simply applying S I SQ PlD control l ors to oach of 
tho n procoss var i ab l os, and 

2. th e r e i s no s e t po i nt in any of a(k-\1),.a(k t 2), . .\, and a(kj) t which can avo i d th e unwant e d sudden 
change in CO when SP changes yyith time,.. 

Tho noxt problom of ■ dosigning tho optimal PID controllor i s to find tho bost va l u es for tho P I D 
param e t e rs K1, K2, . and Kj. An optimization based m e thod for so l v i ng th i s probl e m cons i sts of th o 
fo l lowing four stops: ' 7 -'*; / i;: 

1 . Convort the PID contro l oquat i on into d i scr e t e t i m e form i f i t is not i n d is cr e t e t i m e form. 

2. Bu il d a d i screte time l i n ea r mod el for th e proc e ss th a t is to bo contro ll ed by tho sa i d P I D contro l lor. 

3. Form tho discroto t i mo c l o s od l oop transfer funct i on from sa i d voctor SP to sa i d v e ctor PV. 

A. F i nd tho bost P I D paramotors by using an opt i mization a l gor i thm wh i ch m i n i m i zos tho l argest 
modu l us of a ll po l os of tho discr e t e t i m e c l os e d l oop transfer funct i on obta i ned at stop 3, wh e r e the 
modu l us of a polo i s defined to bo tho absolute va l ue of tho comp l ex number wh i ch represents tho po l o. 
I f tho PID parameters aro subjoct to somo constra i nts, thon a constrainod optim i zat i on a l gorithm can bo 
usod wh i ch min i mizos tho largost modulus of a ll polos of tho d is croto t i mo c l osod l oop tran s for funct i on 
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obta i nod at stop 3 and at tho samo t i mo guarantoos that all usor proscribod constra i nts on tho P I D 
paramotorc aro caticf i od. 

P I D controllers w i th th ei r param e t e rs so obtain e d guarant ee that PV can track SP qu i ck l y 

SPECIFICATION 

TITLE OF THE INVENTION 

Methods of Designing Optimal Linear Controllers 

CROSS-REFERENCE TO RELATED APPLICATIONS 

Domestic priority data: This application is a 371 of PCT/IB01/01002 06/07/2001 
Foreign applications: June 20. 2000 fCAl Canada 2;31 1 ,268 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT 
Not Applicable. t 

INCORPORATlON-BY-REFERfeNCE OF MATERIAL SUBMITTED ON A COMPACT DISC 
Not Applicable 

BACKGROUND OF THE INVENTIONS 
Technical Field of the Invention 

This invention relates to the design of optimal linear controllers and PID controllers- 
Description of Related Art 

Figure 1 shows a process 1 controlled by a PID (proportional-integral-derivative) or a linear 
controller 2. The PID controller or linear controller 2 means a linear system, usually implemented in a 
computer or as an electronic circuit, that receives the process variable signal y(k) and set-point signal 
r(k) as its inputs and calculates the controller output signal u(k) according to a PID control equation or a 
linear control eguation. where k is the integer discrete time variable. For a multiple-input and 
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multiple-output (MIMO) process, the process variable v(k) is an n-dimensional vector with each of its 
components being a scalar process variable and the controller output u(k) is an m-dimensional vector 
variable with each of its components being a scalar controller output, where n and m are positive 
integers. It is desired that the performance of the controller should be such that, after the controller 
output signal u(k) is sent to the process 1. the process variable v(k) should approach the set-point r(k) 
guicklv and smoothly as the time variable k increases, where r(k) is an n-dimensional set-point signal 
(also known as the reference signal or the command signal, etc.) 

There are many types of PIP controllers, depending on the use of different types of PIP control 
eguations. All types of PIP controllers can be viewed as special cases of linear controllers. By definition, 
a discrete-time linear controller means a controller whose output u(k) is a linear function of past 
controller outputs, current and past process variables; and current arid past set-point variables, as 
shown in Figure 2. This definition is well known to any person skilled in the art (seev e.g., Figure 2.4 and 
Assumption 2.3 on page 29 of the cited book "Linear Controller Pesign" authored by Boyd and Barratt 



and published by Prentice Hall in 1991). 




Once the structure of its control equation is properly selected, the performance of a PIP controller 
or a linear controller depends mainly on the choice of coefficients in the control eguation. The 
adjustable numbers in the coefficients are called tuning parameters (also known as tuning coefficients, 
tuning gains or tuning constants' ^tc.) v . . 

How to properly choose tuning parameters for a PIP controller is a problem that has attracted a lot 
of studies ever since PI P Controllers^ ^became widely used in industry in the early 1940s. The 
Ziegler-Nichols tuning methods developed by Ziegler and Nichols in 1942 (see the Ziegler and Nichols 
reference), either in their original form or in some modification, are still widely used in industry. Other 
model-based optimization methods choose the tuning parameters by minimizing some well-known 
control performance index such as the integrated absolute errors (IAE). the integrated sguared errors 
(ISE), the integrated time absolute error (ITAE). etc. (see US patent 5.453,925). However, practice 
shows that all these methods often lead to the undesired oscillatory control results. 

Many tuning methods are based on the continuous-time transfer function analysis (see, for 
example, the cited US patents 6.434.436 B1 to Adamv et al.. 5.866.861 to Raiamani et al.. 4.539.633 to 
Shigemasa et al., 4,563,734 to Mori et al.. the cited research papers by Celentano et al. t Zhuang et al.. 
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Dorf et al.. Katebi et al.. and Saeki et al.. and the cited book "Linear Controller Design" authored by 
Bovd et al.). Therefore these methods cannot guarantee the optimal performance of controllers in the 
discrete-time domain. This is a real problem since nowadays almost all PIP or linear controllers are 
implemented in digital computers and therefore they are working in the discrete-time domain, not the 
continuous-time domain- 
Other tuning methods in discrete-time domain cannot guarantee that the largest absolute value of 
all the poles of the closed-loop transfer function is minimized. Therefore they cannot guarantee the 
optimal performance of the controller (see, e.g., the cited US patent No. 5,680,304 to Wang et al. and 
the cited research paper by Yamamoto et al.. etc.). 

BRIEF SUMMARY OF THE INVENTION 

This invention chooses the tuning parameters in a PI P controller or a linear controller by minimizing 
the maximum of absolute values of all poles of the discrete-time closed-loop transfer function from said 
set-point r(k) to said process variable v(k) subject to, if any, user-specified constraints on one or more 
of the tuning parameters. When the tuning parameters are chosen this way, the PIP controller or linear 
controller can guarantee that the process variable v(k> tracks thS set-point r(k) smoothly and quickly as 
time k increases. %A\. %y 




BRIEF PESCRIPTION bF THE PRAWINGS 

Figure 1 shows a process. 1 controlled by a PIP controller or a linear controller 2. The PIP controller 
or linear controller is in discrete-time foifrn. 

Figure 2 shows a PIP or a linear controller. By definition, the controller output u(k) of a linear 

controller is a linear function of u(k-1). u(k-2) u(k-a), y(k), y(k-1) y(k-b), r(k), r(k-1), r(k-2) 

r(k-c). wherein a>0. b and c are integers (This definition is well known to any person skilled in the art 
and is given in Figure 2.4 and Assumption 2.3 on page 29 of the cited book "Linear Controller Pesion" 
authored by Boyd et al.) 

PETAILEP PESCRIPTION OF THE INVENTION 
From now on it is always assumed that: 
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(3) The open-loop discrete-time transfer function of the process 1 is known, and 

(4) The structure of the control equation of the PIP controller or the linear controller has already been 
properly chosen and is already in discrete-time form. 

It is then easy for anyone skilled in the art to find the discrete-time closed-loop transfer function from 
the set-point r(k) to the process variable y(k). as shown in Figure 1. This invention chooses the best 
values for the tuning parameters in a PIP controller or a linear controller in such a way that the largest 
absolute value of all poles of said discrete-time closed-ioop transfer function from said set-point r(k) to 
said process variable y(k) is minimized subject to. if any, user-specified constraints on one or more of 
the tuning parameters. This choice guarantees that the process^ variable v(k) tracks the set-point r(k) 
smoothly and quickly as time k increases. 

In the Claims: 

Please cancel all claims in record and replace them with the new claims 22-33 (underlined) as follows: 

1. AMIMO (mu l t i p le-i nput mu l t i p le- output) P I D coritroller \^iGh has an n - d i mons i onal process var i ab le 
vector PV with tho n procoss variab l es Pyi, PV2, . :% , and PVh bo i ng i ts first, second, . . . , and n th 
componont, an n-dimonsional sot point vootpr SP with the n sot po i nts SP1, SP2, . . . t and SPn bo i ng 
i ts first, socond, . . . , and n - th comppnont; qnd an m-dimonsional contro l lor output voctor CO w i th tho m 
contro l ler outputs CQ1 ; 002, . . . , and QOm being its first, s e cond, . . . , and m-th component, whoro m 

and n a*e posit i vo 'intogorc, fe; and m wh i ch the P\D — Gontrol e quation is 

CO(k)-CO(k-1)+K1*SP(k)*T^ = . = +Kj*a(k,j), whoro k is tho d i scroto t i mo, T i s tho 

sampl i ng p e r i od, j i s a pos i t i vo i ntogor, K1, K2, . . . , Kj are m by n P I P par a m e t e r s , a(k,1)=[ - PV(k)]*T, 
and a(kj) - [a(kj - 1) - a(k - 1 ,j - 1)]/T for j* or -2. 

2. An M I MO P I P contro ll or of c l a i m 1, i n which tho m by n P I P paramotor s K1, K2, . . . , and Kj aro 
obtain e d by us i ng an opt i mization a l gor i thm wh i ch minim i zes th e largest modu l us of a ll po le s of th o 
d i scroto t i me c l osed l oop tr a nsf e r function from sa i d SP to sa i d PV. 

3. An M I MO P I P contro ll or of claim 2, whoroin tho sa i d opt i m i zat i on a l gorithm is a constra i nod 



PACE 17/32 * RCVD AT 7/12/2004 1:26:18 AM [Eastern Daylight Time] * SVR:USPTO-EFXRF-1/0 * DNIS: 8729306 * CSID: 7804245074 * DURATION (mm-ss): 17-14 



07/11/2004 23:25 



Dan Liu 



PAGE 18/32 



Application/Control Number 1 0/069,364, Art Unit 2121 Page 7 

optim i zat i on a l gor i thm wh i ch m i nimizos tho l argost modu l us of a l l polo s of th e d i scroto t i mo G l osod l oop 
transfor function from s a i d SP to s a i d PV and at tho samo timo guarantee s that tho usor proscribod 
constra i nts on tho P I D parameters aro satisfied 

4 . An M I MO P I D contro ll er of c l aim 1, whoro i n somo or a ll of tho torms K2*a(k,2), K3*a(k,3), , . . , and 
Kj*a(k,j) that appoar on tho right - hand s i de of tho P I D contro l oquat i on aro romovod, for oxamplo, a P I D 
contro llo r with its contro l oquation boing CO(k)-CO(k - 1HK1*SP(k)*T+Kra( - k,1)-CO(k 1) 
+K1*[SP(k) - PV(k)]*T, wh i ch i s a l so ca ll od a 1 - on l y contro ll er, and a P I D contro l ler with its contro l 
oquation — bemg — CO(k)-CO(k - 1)+K1*SF(^^ 1)+K1*[SP(k) - PV(k)]*T K 
2*[PV(k) - PV(k - 1)], wh i ch i s a l so cal l od a PI control l or, otc. - t 

5. An M I MQ P I D contro lle r of c l aim 4 , wh e r ei n th e remaining PID parameters are obtained by us i ng an 
opt i mizat i on a l gor i thm which m i n i m i zos tho l argost modu l us of a ll po l os of tho d i scroto t i mo c l osod l oop 
transfor funct i on from sa i d SP to sa i d PV. 

6. An M I MO P I D controllor of iclaim S-jyvheroin the sa i d opt i m i zat i on a l gorithm i s a constra i ned 
opt i mizat i on a l gor i thm wh i ch m i n i mi(^,thQ > ^r-gpst modulus of a l l po l os of tho d i scroto t i mo c l osod l oop 
transfor funct i on from said SP to ggid PV and at tho samo tim e guarantoos that the u s er proscr i bod 
constra i nts on th e PID parameters are Satisfied. 

7. A P I D contro ll or of claim 1, vVhorpin said PV, sa i d SP, sa i d CO, and said P I D paramotors aro alt 
sca l ars, and m~n~1. 

8. A P I D contro ll or of claim 2, whoroin sa i d PV, said SP, sa i d CO, and sa i d P I D paramotors aro a lt 
sca l ars, and m-n-1. 

9. A P I D contro ll er of c l a i m 3, whoroin sa i d PV, sa i d SP, said CO, and sa i d P I D paramotors aro a ll 
scalars, and m-n-1 . 
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10. A F I D contro ll er of claim 4, whoroin sa i d FV, sa i d SF, said CO, and s aid P I D paramotors aro a l l 
s ca l ars, and rrrn-1. 

11. A P I D control l er of cla i m 5, wher e in s ai d PV, sa i d SP, sa i d CO , a nd sa i d PIP par a m e t e rs ar e a l l 
scalars, and m-n~1 . 

12. A F I D contro ll or of c l aim 6, whoro i n sa i d FV, sa i d SF, sa i d CO, and sa i d P I D paramotors aro a ll 
scalars, and nrrn-1. 

13. Amothod of f i nd i ng tho opt i ma l P I D paramotor s or tuning constants for any P I D contro ll or or a li noar 
contro lle r by m i n i miz i ng tho largost modu l u s of a ll pol e s of the d i scr e t e- tim e c l os e d - loop transfor 
function from s e t po i nt SP to process var i ab l e PV. 

14. A mothod of C l aim 13, whoroin tho optimal PID ;paramc^ors pr ; tun i ng constants aro obtainod by 
m i n i m i zing th e l arg e st modu l us of the discrete-time clos^d^lpop transfer funct i on from said SP to sa i d 
PV s ubject to tho constraint that spmo or ail of tho PiD paramotors or tun i ng constants aro w i thin usor 
spoc i fiod admiss i b l o rangos. ^ $ ^ 

15. A PID contro ll er of Claim 1 or any FID controller or li n e ar control ler, wh e r ei n th e PID p a ramotors or 
tun i ng constants aro obtaindcl such that fc: 

(1) tho max i mum roa l part of a ll polos of tho cont i nuou s- time c l osod -l oop tran s for funct i on from said SP 
to sa i d PV i s m i nimiz e d, poss i bly subjoct to any usor sp e c i fi e d con s traints on tho P I D param e t e rs or 
tuning constants, such as som e or al l of eleme nts i n th e P I D paramet e rs or tun i ng constants shou l d 
bo w i th i n usor spoc i fied rangos, whoro tho sa i d P I D paramotors or tun i ng constants can bo matr i ces 
or sca l ars, or 

(2) the- max i mum magn i tude of som e or a l l ele m e nts i n the PID paramet e rs or tun i ng constants i s - 
m i n i mized, subjoct to tho constra i nt that tho max i mum r o a l part of a ll po l os of tho sa i d 
cont i nuou s- t i mo c l osod -l oop transfor funct i on from s a i d SF to sa i d FV is not l argor than a usor 
spoc i fiod numbor. 
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22 A method for determining the optimal tuning parameters in a linear controller wherein 

a) said controller receives an n-dimensional process variable signal v(k) from a process and an 
n-dimensional set-point signal r(k). calculates an m-dimensional controller output u(k) 
according to a linear control equation, and sends said u(IO to said process, where k is the 
integer discrete time variable and n and m are positive integers; 

b) said tuning parameters are the adjustable numbers in the coefficients in said linear control 
equation that are to be determined; and 

c) said method finds the optimal values for said tuning parameters by minimizing the maximum of 
absolute values of all poles of the discrete-time closed^oop transfer function from said 
set-point r(k) to said process variable v(k). 

23 A method as in Claim 22, wherein said minimization of the maximum of absolute values of all poles 
of said discrete-time closed-loop transfer function is subject to user-specified constraints placed on 
one or more of said tuning parameters. 

24 A method as in Claim 22. wherein said controller output u(k) ^ u(k-1) ± K i *r(k)*T ± K i *a(k.1) + 
K /a(k.2) + ... + K P *a(k.p). wherein k is the discrete time variable. * is the multiplication operator T 
is the sampling period, p is a positive integer, K i , K 2 , and K p are m by n tuning parameters. a(k.1) 
= \-v(k)VT. and a(k- p) = ra(kfel) - a(k-1.p-1)l 7 t for p> or =2. 

25 A method as in Claim 23, wherein said dontroller output ufk) = u(k-1) + K i »r(k)*T + K i *a(k,1) + 
K?*a(k.2) + ... +K n *a(k, p> wherein k the discrete time variable. * is the multiplication operator. T is 

the sampling period, p is a pd^itive iriteger. K i , K 2 and K c are m by n tuning parameters. a(k.1) = 

hv(k)rT. and a(k. p) = fa(k.p-1) - a(k-1 .p-1)l / T for p> or =2. 

26 A method as in Claim 22, wherein said linear controller is a PIP (proportional-integral- derivative) 
controller. 

27 A method as in Claim 23. wherein said linear controller is a PIP controller. 

28 A linear controller as in Claim 22 with its tuning parameters determined therein. 

29 A linear controller as in Claim 23 with its tuning parameters determined therein. 

30 A linear controller as in Claim 24 with its tuning parameters determined therein. 

31 A linear controller as in Claim 25 with its tuning parameters determined therein. 
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32 A PIP controller as in Claim 26 with its tuning parameters determined therein. 

33 A PIP controller as in Claim 27 with its tuning parameters determined therein. 

In the Abstract 

Please cancel the old Abstract in record and replace it with the new Abstract (underlined) as follows: 

Mothods of dosign i ng tho structuro of multip l o -i nput mu l tip l o - output (M I MO) P I P contro l lors and 
methods of f i nding th e opt i ma l va l ues for th e M I MQ P I P param e t e rs ar e d i sc l os e d. Tho opt i ma l va l u e s 
of MIMO P I P param e ter s aro obta i nod by us i ng an optim i zat i on algor i thm which m i n i mizos tho l argoct 
modu l us of all po l os or tho d i scroto timo closod l oop transfer function from sot po i nt SP to propose 
var i abl e PV, with or w i thout usor proscr i bod constraints on th e PIP paramotors. Mothods of dosign i ng 
th e structur e of s i ngl e -input s i ngl e- output (S I SO) PIP cqntrolbrs and methods of finding t il o opt i ma l 
valuos for S I SQ PIP paramotors aro also disclosed as spocial MIMO P I P contro ll or casos. 

Methods of designing optimal discrete^time PIP (proportiona controllers and linear 

controllers are disclosed. The optimal values of the tuning parameters in a PIP controller or a linear 
controller are determined by miriimizing the maximum of absolute values of all poles of the 
discrete-time closed-Jbop transfefAunctidti from thb set-point to the process variable subject to. if any, 
user-specified constraints on one or more of th^ tuninq parameters. 

In the drawings: 

Please add the attached Figure 1 arid Figure 2. 
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